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1. Introduction 



Two-dimensional spin-^ non-relativistic quantum systems with magnetic fields 
are described by the Pauli operator. For regular magnetic fields the Pauli operator 
is usually defined as the square of the Dirac operator. However, for more singular 
^C ' magnetic fields, such as the delta field, an Aharonov-Bohm (AB) solenoid, gener- 

ates (see |2]), the situation is more complicated. Then there are many self-adjoint 
^X ■ extensions of both the Dirac and the Pauli operator, originally defined on smooth 

functions with compact support not touching the singular points. The different 
^sO ■ extensions describe different physics, and it is not clear which extensions describe 

the real physical situation. 

We consider the case of a magnetic field consisting of finitely many AB solenoids 
and a smooth field with compact support. Up to now only two different Pauli 
extensions have been studied for this type of magnetic field (see ^1 E]), both 
defined via a quadratic form. Since the Pauli operator classically is the square of 
the Dirac operator it is motivated to study the self-adjoint Pauli extensions that 
can be obtained in this way. This is the purpose of the present paper. 

A natural property to expect from the Pauli operator is that it transforms in an 
(anti)-unitarily way when the sign of the magnetic field is changed and the spin-up 
and spin-down components are switched. This property is usually called spin-flip 
invariance, and we want to answer the question of which Pauli operators defined in 
different ways satisfy it. 

Another natural property to expect is the possibility to approximate our operator 
with operators corresponding to regular magnetic fields. In |Hj it was studied which 
Pauli extensions are possible to approximate, and the condition was expressed in 
terms of the asymptotics of the functions in the domain at the singular points. 

One of the ways to control the properties of self-adjoint extensions consists in the 
analysis of the dimension of the kernel of the operator. The dimension of the kernel 
of the Pauli operator is usually given by the Aharonov-Casher formula (see @|). 
The formula has been proved in different settings, see |51 1131 IT7] . Recently this 
formula was also proved for one of the extensions for a very singular magnetic field 
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(containing the case with AB solenoids) in ^5| • Another extension was introduced 
in ^21, and in ^B] an Aharonov-Casher type formula was established for that 
extension. 

The Dirac operator with strongly singular magnetic field has been studied before 
in 1^1 [7| E3 1151 I2(J| . In JS] a formula for the dimension of the kernel of the 
Dirac operator was proved for two different asymmetric self-adjoint extensions (i.e. 
those with different behavior of spin- up and spin-down components), and it was 
demonstrated that, in fact, this dimension may differ for self-adjoint realizations, 
although, each of them seems to be quite natural. These extensions are closely 
related to the ones introduced in 6 . In both these articles the magnetic field is 
the same as the one we consider (the one in |S] does not have the regular part), 
with the addition of even more singular terms containing derivatives of the delta 
distributions. Using gauge transformations one can dispose of these derivatives. 

In |2()j it was proved that one of these asymmetric Dirac extensions can be 
approximated by operators corresponding to more regular magnetic fields. However, 
this extension, as we show, lacks the property of being spin-flip invariant. 

In Section [2] we study the Dirac operator. In order to be able to treat the 
general case, we need first to repeat in details the description of all self-adjoint 
extensions corresponding to only one AB solenoid, given in |20j . To extend this 
description to the case of several solenoids, we use the glue-together procedure, 
proposed in 2 . This procedure enables one to glue together operators with different 
behavior at different points; however it is natural to restrict oneself to the case of 
the same behavior at all points. After that we check which extensions are spin-flip 
invariant and finally we prove a formula for the dimension of the kernel of the Dirac 
extensions. 

In Section [3] we consider the Pauli operators that are the square of some self- 
adjoint Dirac operator defined in Section|21 We show exactly which Pauli extensions 
are obtained in this way, in terms of the asymptotics of functions in the domain 
of the Pauli operator at the points where the singular AB solenoids are located. 
We also find an Aharonov-Casher type formula for these Pauli operators. It turns 
out that there are only two of them that have zero-modes. These two extensions 
are very asymmetric though, admitting singularities in one component only, which 
looks rather non-physical. All the other extensions have singularities in both the 
spin-up and spin-down components, and they are coupled. 

It turns out that the Pauli operator studied in |12j is a mixture of these two asym- 
metric extensions, admitting different physical situations at different AB solenoids. 
This has to do with the normalization of the AB intensities. In that article the AB 
intensities were chosen to be normalized to the interval [—1/2, 1/2). If any of the 
intervals (0, 1) or (—1,0) would have been chosen instead, then the Pauli operator 
would have been one of the asymmetric ones studied in this article. Both these 
asymmetric operators have the advantage that they describe the same physics at 
all AB solenoids. In the end of the article we present a discussion of the properties 
of the self-adjoint Pauli extensions with respect to different ways of normalization 
of AB intensities. 



2. The Dirac operator with singular magnetic field 

The goal in this section is to describe the self-adjoint Dirac operators corre- 
sponding to a magnetic field consisting of several (but finitely many) AB solenoids 
together with a smooth field, and to find an Aharonov-Casher type formula for 
the dimension of the kernel of these self-adjoint operators. Let us introduce some 
notations that will be used throughout the article. We identify a point (£1,3:2) in 
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K with the complex number z = X\ +1x2, and we will often write z in polar coordi- 
nates, z — re l9 . Sometimes it will be convenient to use the polar coordinates rje l6j 
with Zj as the origin. The magnetic field will consist of a regular part Bo <G Cq (M 2 ) 
and a singular part consisting of n AB solenoids located at the points A = {zj}™, 
so that the magnetic field B has the form 

n 

B(z) = B {z) + J2 2 ™j5 Zj . (2.1) 

i=i 
Owing to gauge equivalence (see [20]) we can assume that all the AB intensities 
o.j (fluxes divided by 2-k) belong to the interval (0,1). All derivatives will be 
considered in the distribution space 2?'(C \ A). We will denote by h a magnetic 
scalar potential satisfying Ah = B. The magnetic scalar potential is uniquely 
defined modulo addition of a harmonic function. We will use the scalar potential 
-if n n 

h(z) = — B (Olog\z-C\dX(0 + ^2a j log\z-z j \=h (z) + ^2h j (z), (2.2) 

where dX is the Lebesgue measure. The actions g±, which will be used to describe 
how the Dirac operator acts, are defined by 



q%u = -2i£ (ue- n ) e" 



and 



q h u = -2i-§- z (ue h ) e~ h . 

These actions q+ and q^ are usually called the spin-up and spin-down actions, 
respectively. The Dirac action is given by 



h _ I q'l 



To be able to describe the self-adjoint Dirac operators with several AB solenoids we 
first study the self-adjoint extensions of the Dirac operator with one AB solenoid, 
originally defined on smooth functions with compact support not touching the 
singular point. 

2.1. The Dirac operator with one AB solenoid. The case of one AB solenoid 
has been studied before (see |1()U2()| ). and we just sketch the way it was done since 
we need the detailed information about these extensions for our further analysis. 
We let the AB solenoid have intensity a = ot\ £ (0, 1) and be located at the origin. 
We will describe all self-adjoint extensions of the Dirac operator originally defined 
on (C^M 2 \ {0})) 2 . Let us drop the superscript h of q± and d h in this subsection. 
The actions q± can in this case be written as 

9 + = "Hal " 2l)= - e {8-r + rm--) ^ 

and 

„. / d a\ iff f 9 id a\ . ,. 

The minimal Dirac operator S m in, obviously symmetric, is defined by 
2)(D m i„) = (C °°(IR 2 \{0})) 2 , 
£W/> = Dip, for ip G D(S m ; n ). 

It can be seen that 2) m in has deficiency index (1,1), and the deficiency spaces 
X± = ker(!D min =F i) are spanned by 

«*<•*") -(^rlD- (2 - 5) 
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Denote by U any unitary operator from x+ to \—- Then U takes £+ to e jr £_ for 
some r € [0, 27r). According to the theorem of Krein and von Neumann, described 
in |5], all self-adjoint extensions of !D m j n can be parameterized by r as 

D(S) r ) - {iP e {L 2 (R 2 )) 2 : if, = -00 + M (£+ + e lr £_),0 o € D(S m in), /* € C} 

(2.6) 

® T ip = OVo + «M<£f - e ir £-), for V € , D(S) T ). (2.7) 

It is also possible to describe the self-adjoint extensions by studying the asymptotic 
behavior of the functions in the domain at the origin. To see this, let us define the 
linear functionals c_ a and c^_ 1 on T>(T) T ) as 

c ± M±) = limr"— / ij}±{re ie )d6, and (2.8) 

r— 2-7T J 

ctxW-i) = limr 1 -"^ f \ ± (re l6 y e d6. (2.9) 

For ip = ip + n(£ + + e""£_) in 1>(D T ), where Vo € D(S min ), applying _these 
functionals gives no contribution from ?/>o since the limit of functions in D(Smin) 
tends to zero at the origin. Let us introduce the notation <r{a) — Y(a)2 a . Using 
the asymptotics for the Bessel functions we get 

c±M+) = o 

cl a ty_) = £(l-e*>(a) 

a ( 2 - 10 ) 

4_ 1 (V + ) = f(l + e l >(l-a) 

c"_ 1 (V-) = 

for such functions -0 € ©(S 1 "). Here /i is the same constant as in (|2.6|) . An 
equivalent description of all self-adjoint Dirac extensions is 

£(S T ) = {V; G (L 2 (M 2 )) 2 : ^ e (L 2 (M 2 )) 2 ; (2.11) 

c «-iO+) . ,, /oN o-(l-q) 

_ , . : = icot(r/2 --— , and 

C-M-) a \ a ) 

c±«W>+) = c-.iftM = o}, 

D T ilj = <)ip, for V> S D(S) t ). (2.12) 

2.2. The Dirac operator with several AB solenoids together with a reg- 
ular field. In this subsection we are going to study the Dirac operator for a mag- 
netic field consisting of a finite number of AB solenoids together with a regular 
background field. We will use the same method as in to glue together the dif- 
ferent self-adjoint Dirac operators corresponding to only one AB solenoid and the 
self-adjoint Dirac operator corresponding to the regular magnetic field. 

Note here that we do not study all self-adjoint extensions but only the ones that 
are subject to the natural locality principle; which says that the asymptotics at a 
singular point of functions in the domain should be independent of the asymptotics 
at the other singular points. 

We start by defining the Dirac operator with two AB solenoids together with a 
smooth field. The general case does not give any extra difficulties. Let the magnetic 
field B consist of a smooth field Bq with compact support and two AB solenoids 
located at z\ and zi with intensities ol\ and a 2 , 

B{z) = B (z) + 27rai<5 Zl + 2Tra 2 5 Z2 . (2.13) 
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In this case our scalar potential h can be written as 

h(z) = ho(z) + hi(z) + h 2 (z) = — (log | ■\*Bo)(z) + ailog\z-zi\ + a 2 log\z-Z2\. 

Zn 

From the previous section we have self-adjoint Dirac operators £> Tl '' 11 and J) T2 ''' 2 
corresponding to each of the AB solenoids separately. Let us drop the parameters 
Ti and T2 from the superscripts. So, for example, when we write 1) hl we mean some 
self-adjoint extension with one AB solenoid located at z\. 

Let (fij G C^°(R 2 ), j — 1,2, be equal to 1 in a neighborhood of Zj and have 
small support not touching a neighborhood of Zk, k ^ j and < ipj < 1. Let 
tpo = I — ipi — (fi2- We denote by Ejk the set suppi^,- n suppipk- 

Let us introduce the multiplication operators V hj as 

vhj - 2 iA- ~t j 

\dz'h u 

Note that V h ° is bounded from (L2) 2 to (L2) 2 . For j ^ we will be sure to apply 
the operators V hj only on functions being zero in a neighborhood of the singular 
points Zj. 

Definition 2.1. The Dirac operator T) h corresponding to the magnetic field B 
in H2.13|l is defined as 

V(D h ) = {V G (i 2 (K 2 )) 2 : <Pjil> G B(S' l O, j = 0,1,2} 

and 

I)V =(D' 10 + V hl + V h2 ){Lp Q ^) 

+ (T) hl +V ho +V h2 )(tpnlj) 
+ (® h2 +V ho +V hl )(ip 2 i;) 
for ^ G 2)(D''). 

It is easily verified that the definition is independent of the partition of unity 

1 = ifio + (fii +lfi 2 - 

Theorem 2.2. The Dirac operator 'D h is self-adjoint. 

For the proof, we need some lemmas. 

Lemma 2.3. The Dirac operator T) : (D2) 2 — > (L2) 2 without any magnetic field is 
a self-adjoint operator with the Sobolev space (H 1 ) 2 as domain. 

Proof. See EI]- □ 

Lemma 2.4. The Dirac operator 2) ° corresponding to the magnetic field Bq is 
self-adjoint with domain (H 1 ) 2 . 

Proof. The operator T) h " can be written as 1) h " = T) + V h " and the multiplication 
operator V h ° is relatively bounded with respect to J) with relative bound zero, so 
the lemma follows from the Kato-Rellich theorem. □ 

Lemma 2.5. Let T be a bounded operator from (L2) 2 to (H 1 ) 2 and let V be a 
function, V(z) — > as \z\ — > 00. Then the composition VT is compact in (L2) 2 ■ 

Proof. For n = 1, 2, ... we write V as V = V n + V n , where 

4) \0 W(z)\<±. 

The functions V n all have compact support, so the operators V n T are compact. But 

\\V n T -VT\\ < i||T|| for all n = 1,2,..., so VT is also compact. D 
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Remark. Lemma 12.51 is also true for 2 x 2 matrix valued functions V where all 
components tend to zero at infinity. It also holds if T is bounded from L 2 to H 1 . 

Lemma 2.6. Let / s e R and let ip £ C^ D (M 2 ) with zero in its support. Then 
the operator ipR is compact, where R — (2) T + is) -1 and T) r is any self-adjoint 
extension of the Dirac operator corresponding to one AB solenoid (which is assumed 
to be located at the origin). 

Proof. First, ipR is compact if and only if ipR((pR)* = ipRR*ip is compact. To 
show that <pRR*tp is compact, it is sufficient to show that (pRR* is compact. 

The operator RR* is equal to ((S T ) 2 + s 2 )~\ Note that (£> r ) 2 is a self-adjoint 
Pauli operator corresponding to the same magnetic field (see Section |3| for a dis- 
cussion of the Pauli operators that are the square of some Dirac operator) . If we 
denote by ty any other self-adjoint Pauli operator corresponding to this magnetic 
field, then by the Krein resolvent formula (see 0) the resolvents of (J) 1 ") 2 and ty 
differ by a finite rank operator. Thus, it is enough to show that <p(?p + s 2 ) -1 is 
compact for a convenient choice of self-adjoint Pauli extension <$. Let us choose 
*P to be the Friedrichs extension. The functions in the domain of this extension ^3 
vanish at the origin so 



Hj 

where H is the Friedrichs extension of the Schrbdinger operator corresponding to 
the same magnetic field (see ^2 f° r a discussion of this). Hence it is enough to 
show that (p(H + s 2 ) -1 is compact. 

Let Hq = —A be the Schrodinger operator corresponding to no magnetic field. 
Then, by the diamagnetic inequality (see ^B]) it follows that \tp(H + s 2 ) _1 u| < 
</?(i/o + s 2 ) _1 |u| (pointwise) for all u £ L 2 . This inequality implies that ip(H + s 2 ) -1 
is compact if (p(H + s 2 ) -1 is compact (see [TTUTT?] ). 

The compactness of ip(Ho + s 2 )^ 1 follows from Lemma \'2. 51 since (Hq + s 2 ) -1 is 
bounded from L 2 to H 1 . □ 

Lemma 2.7. The operator 2) is symmetric. 

Proof. Let ip and tp belong to T>(S) h ). Then 

(33 V, i>) = /(2» V) ' i>d\(z) = f((cp + Vl + pa)3>V) ■ 4>d\(z). 
The symmetry now follows from integration by parts, noticing that (for example) 

<Po(® h ° + V hl + V h2 ){<p <l>) ■ HX(z) = f ipotp- {T) h <> + V h * + V^){ipo4>)d\{z) 
and 



^i(2T° + V 11 + V n2 )((p tp) ■ ^dX(z) = I (S)' 10 + V 111 + V**){ywl>) ■ ipi^dX(z) 



= / <p ip ■ (£> ftl + V h ° + V h ^)(ip 1 Tp)d\(z) 

and similar for the other terms. Adding all terms we see that T) h is symmetric. □ 

In the following lemma we look at our operator as acting from its domain T){T) h ) 
considered as a Hilbert space equipped with graph norm 

2 h = IKS)' 10 + v hl + v h2 ){ Vo m 2 + IK®' 11 + v ho + T/^X^VOH 2 



is 

+ \\(D h2 +V h ° +V hl )(^ 2 ^)\\ 2 + 
to (L 2 ) 2 . 
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Lemma 2.8. Let £ a £ K be fixed. The operator 5) h + is : (£>(£''), || • || S h) -» 
(L2(R 2 )) 2 is a bounded Fredholm operator with index zero. 

Proof. First, it is clear that S" + is is bounded from the domain space with 
graph norm. To show that 1) h + is is a Fredholm operator, it is enough to find a 
left and a right parametrix (see pQ). We start by finding a right parametrix. Let 
Rj denote the resolvent Rj = { r £) hj + is) -1 , j = 0,1,2, and define the operator 

R : (L 2 f -► (L 2 ) 2 as 

Ru = LfoRou + ipiRiu + ip2R 2 u 1 for u £ (L2) 2 . 

For u £ (L2) 2 we have ipjifkRjU £ (H 1 ) 2 and being zero in a neighborhood of the 
singular point(s) if j =/= k. Thus 

(j)fc* + V^frjtpkRju) = (S^ + V h *){<p m R jU ), j £ k. 

From this it follows that, if u £ {L2) 2 and j, k, I £ {0, 1,2, j ^ k =^ /}, then 

(£' 1 + is^Rju) = <p 3 u + {{V hk + V hl )<pj + ®{ip 3 ))R.jU. 

Thus 

(D h + is)Ru = u + K R u 
where Kr : (L2) 2 — ► (L2) 2 is the operator 

K R u = ((V hl + V h2 )<p Q + D(ip )) R Q u 
+ ((V ho + V' 12 )^ + D(ifi)) Rm 
+ ((V ho + V hl ) V2 + ®(<pa)) R 2U- 

Kr is compact. Indeed, the first term is compact according to Lemma l2.5l since Rq is 
bounded from (L 2 ) 2 to (if 1 ) 2 and the matrix- valued function (V hl +V h2 )(p Q +'B(ipo) 
tends to zero at infinity. The other two terms are compact by Lemma J2.6I Hence 
Kr is compact, so R is a right parametrix. 

The operator L — Rq<po + Ri <fii+R2<fi2 is a left parametrix. Indeed, a calculation 
similar to the one above shows that 

L(D h + is)if) = V + K L iP 

where 

K L ij, = R {(V h ^ + V h2 )cp - JDftx,)) V 

+ R 1 ((V ho +V h2 )i Pl ~®{vi))^ 

+ R2((V ho +V h i)^2-V(¥2))^. 
The compactness of K* L follows in the same way as the compactness of Kr. Hence 
Kl is compact. We see that any of R and L works as a parametrix and hence 
T) h + is is a Fredholm operator. 

To see that 1) h + is has index zero, we note that since 3} with domain (H 1 ) 2 is 
self- adjoint, it has index zero and R s := (2) + is)~ 1 is a parametrix for iD + is. The 
operator 

R- R s = poRo + fiRi + V2R2 - Rs 

is compact. To see this, write R — R s as 

R - R s = {fo - l)iZo + <Pi&i + V2R2 + (Ro - R s )- 

The first term is compact according to Lemma 12.51 the second and third according 
to Lemma 12.61 For the last term we note that Rq — R s = —RoV h °R s . The 
compactness of V h °R s follows from Lemma [2.51 Composition with the bounded 
operator R preserves compactness. Thus R — R s is compact. It follows that 

ind(.R) = ind(i? s ). 
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Since R and R s are parametrices for 1) h + is and £> + is respectively, it holds that 
ind(.R) + ind(S)' 1 + is) = md(R(5) h + is)) = 0, 



and 
Hence 



ind(i? s ) + ind(£ + is) = ind(i? s (D + is)) = 0. 
ind(2)' 1 + is) = - ind(-R) = - md(R s ) = ind(£ + is) = 0. 



□ 



Proof (of Theorem 12 .2JI . We know from Lemma \'2 . 71 that 2) is symmetric, so for 
7^ s £ R we have 

||(D' 1 +zs)V|| 2 = ||S)V|| 2 + s 2 ||V-|| 2 > s 2 "-'-" 2 



It follows that dimker(S)'' + is) — 0. From Lemma f2. 81 we have that T) h + is has 
index zero, so it follows that dimker((D /l )* — is) — 0. Choosing s positive and 
negative respectively gives that the deficiency indices for 1) h is (0,0), so S^ is 
self-adjoint. □ 

2.3. Spin flip invariance. Since the particle we are studying moves only in a 
plane, and the magnetic field is orthogonal to this plane, physically it should be 
no difference if the sign of the magnetic field is changed. This transformation 
has to come together with a flip of the spin-up and spin-down components and a 
normalization of the AB intensities. We say that a self-adjoint extension is spin flip 
invariant if, after applying these transformations, we end up with a (anti)-unitarily 
equivalent operator. We will show that there are only two values of r for which 
the Dirac operator J) 1 "'' 1 is spin flip invariant. Let t — {t\,.. . ,t„) and denote the 
Dirac operator by 1)' r ' h . We will use the linear functionals 

i r 27r 

1 ■ ' J0i 



ct a M±)= limr f— / iP±(r je ^)d0 3l and (2.14) 

4_i(^±) = Kmr 1 -**!- f \ ± (r 3 e^)e^de 3 . (2.15) 

We define anti-unitarily operator V : (L2(M. 2 )) 2 —> (L2(M 2 )) 2 as the spin-flip oper- 
ator that takes (i/>+,'0-) t to (V'-,i/'+)*- 

Proposition 2.9. The operators I) 1 ' '~ h and '£> Th are anti-unitarily equivalent via 
the operator V if and only if for all j = 1, . . . , n we have t'j+Tj = tt or t'^ + tj = 3n. 

Proof. Let j3j = 1 — (Xj be the normalized AB intensities for the magnetic field 
— B that corresponds to D 1 " > . A function ip in the domain of £> T ' '~ h has the 
asymptotics 

H /(I + e^)a(l Mrj'- 1 + O^A 

as z — > Zj for some constant fij £ C. We see that Vtp has the asymptotics 
(1 - e- iT 'i)a{l - aj)rf~ x e- i6 i + 0(r]~ a ')\ 
(l + e-^)a(a 3 )rj a > +0(r^) ) 

Applying the functionals l|2.14|l and (|2.15|) we see that Vtp satisfies 




c 



+ 



-i((^)+) ., , , /n Mi 



itan(r'/2)- 



,((V^)_) VJ/ <t{oj) 



PAULI OPERATOR WITH SEVERAL AB SOLENOIDS 9 

and c~._ 1 ((V'0-)-) = C-a • (O' r, 0+)+) = 0, so the requirements that the domain 
change properly is that 

tan(rj/2) = cot(rj/2), for j = 1, . . . , n. 

We see that Tj/2 and tt/2 — rj/2 must differ by a integer multiple of it. Both tj 
and r' belong to the interval [0, 2tt), so the only possibilities are r.' + tj = tt or 
., = 3tt. D 



t j + 7 - 

Corollary 2.10. TTie operators T) T, ~ h and T) T,h are anti-unitarily equivalent via 
the operator V if and only if for all j = 1, . . . ,n we have tj — n/2 or Tj = 3n/2. 

Proof. Take r' = tj in the previous Proposition. D 

If we let W : (L 2 (1R 2 )) 2 -> (L 2 (R 2 )) 2 be the operator that takes (i/> + ,ip_) t to 
(if)—, ip+y we get some other symmetries if we compose it with the gauge transform 
that only act on the spin- up component. 

Proposition 2.11. The operators T) T '~ h and r D T ' h are unitarily equivalent via 
the operator W composed with a gauge multiplication of e~ 2t ^ j of the spin-up 
component if and only if |r' — tj\ = tt for all j = 1, . . . , n. 

Proof. The proof goes on as in the proof of Proposition ^ . 91 This time the require- 
ment on Tj and t' becomes 

— tan(rj/2) = cot(rj/2), for j = 1, . . . , n 

which gives |r' — Tj\ — it for all j = 1, . . . , n. □ 

2.4. Zero-modes. Let us calculate the dimension of the kernel of 1) h under the 
assumption that Tj = r for all j = 1, ...,n, which means that we assume that 
we have the same physical conditions of the behavior of the particle close to all 
solenoids. Denote by $ the total flux of B divided by 2-7T, that is 

* = T [ B ^ dX ^ = Y I B Q {z)d\(z) +J2°<J- 

J 7T J ^ 

As usual, the definition of the total flux is a matter of agreement, due to the 
arbitrariness in the choice of normalization for AB intensities. The asymptotics of 
e h at infinity and at the singular points A are given by 

h I |z|* \z\ — * 00 , 

[Iz-ZjH z^Zj. 

Remember also that the functions in the domain of T) h must satisfy 

f^.^.fflf^a!, 3 = 1 ,...,„. (,.!„ 

Let {x} denote the lower integer part, that is 

[x\ x > 1 and x ^ N 
{x} = ■{ x — 1 x > 1 and x € N 
otherwise. 

Theorem 2.12. If Tj = t, j = 1, . . . , n then the dimension of the kernel of 1) h is 
given by 

f{|n-$|} ifr = 0, 
dimkerl)' 1 = <^ {|$|} J/t = tt, 

I otherwise. 
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The proof follows the same idea as the original proof by Aharonov-Casher with 
the same changes as in |18| and using the fact that the spin-up and spin-down 
components are coupled if r ^ {0, ir}. 

Proof. We start by calculating the zero-modes as if the spin-up and spin-down 
components were not coupled; so these components are studied separately. 

Let us start with the spin-up component, that is, we consider the solutions to 
q+-0+ = 0- This is equivalent to ■^(e~ h "ip + ) — 0, and thus the function / + = e~ h ip + 
must be analytic in C \ A. The behavior of /+ at the singular points A is different 
for different values of the parameter r, but a pole of order at most {— $} — 1 at 
infinity is allowed independently of the value of t. 

Case I, r = 7r: For ip + to belong to L2, as we see from (|2.16|) . the function / + 
is not allowed to have any poles at the singular points A. Thus, if r = 7r then / + 
may be a polynomial of order at most {— $} — 1. There are as many as {— <&} many 
linearly independent such polynomials. 

Case II, t ^ tt: From (|2.17() we see that a pole of order at most one is allowed 
at each Zj £ A. The calculation in JHj then yields that the dimension is {n — $}. 

Let us now turn to the spin-down component. We look for solutions to the 
equation q~ip- = 0, which is equivalent to finding solutions to ■§^{& h '4>-) = 0. If 
we now let /_ = e ip-, then /_ must be anti-analytic in C \ A, and from the 
asymptotics (|2.16ll we see that /_ may have a polynomial part of degree at most 
{$} — 1 independent of the value of the parameter r. Again we get two different 
cases for the behavior of the functions at the singular points A. 

Case I, r = 0: In this case we see from (I2.17[) that no singular parts for ip_ are 
allowed at A, and hence /_ must have a zero of order at least 1 at each point in 
A. That is we have a polynomial in z of degree {$} — 1 with n predicted zeroes. 
There are {<£> — n} linearly independent polynomials of this type. 

Case II, t^0: Now /_ must be a polynomial in z of degree at most {$} — 1, 
but without any forced zeroes. Thus the dimension of the kernel is {$}■ 

Since the spin-up and spin-down components are not coupled in the cases r = 
and t = 7r the calculations above yield 

dimkerST = { )\ n IJ 

1{|$|} ifr = 7T. 

Let us now assume that r $• {0, 7r}. We should evaluate how the spin-up zero-modes 
match the spin-down zero-modes to satisfy the conditions at the singularities. First 
we note that to be able to have zero-modes both {n — <&} and {<&} must be positive. 
From the calculations in the last two paragraphs of the proof of Theorem 3.3 in JS] 
it follows that /+ must be of the form 

n 

/+(*) = E— 

*—' Z — Zj 
3=1 3 

where r/j 6 C satisfy 

n 

E?7i4=0, for k = 0,l,...,n- {«-$}- 1 (2.18) 

3=1 

and f-(z) must be a polynomial in z of degree at most {$} — 1. Actually, we will 
show that even if the degree of the polynomial /_ is {$} or in some cases {$} + 1, 
all coefficients of the polynomial must be zero. Let us define the natural number 
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in as m = n — {n — <&} — 1 and note that m — L$J ■ Let 

m 

f-( z ) = ^ajZ J . 



From the asymptotics Ij2.17|l we see that 



Vj 



f-(Zj) 



= e- 2h "^ J] (\ Zj - z,r 2tt! ) »cot(r/2) 



¥j 



cr(l — Oj) 



From the requirements (|2.18|l of the coefficients rjj we get 



3 = 1, 



(2.19) 



o = yj^j^ = *cot 



(r/2)£ &,/_&)**, fe = 0,l,...,i 



where 






(2.20) 



We introduce the vector a = (ao ,•■■,< 
Let us also introduce the matrix V as 



V 



/I 1 



V 



1 z 2 



i)' where a, are the coemcients in H2.19J) . 

5 



and the diagonal matrix B having the positive number bj at the j'th diagonal 
position. Then i|2.20[l can be written as 

icot(r/2)VBV*a = 0. 

The matrix VBV* is clearly Hermitian and since B is positive, we can write VBV* 
as (V 'y/~B)(V V~B)* . Hence the null space of VBV* is the same as that of the matrix 
(Vy/~B~)* = \fBV* . Since V* is (a part of) a Vandermonde matrix it has full rank, 
so the dimension of the null space of \f~BV* is zero. Hence the polynomial /_, 
and thus also ?/>_, must be zero. Since the spin-up and spin-down components are 
coupled, it follows that ip+ is also zero. Consequently, dimkerS)' 1 = 0, and the 
proof is complete. □ 



3. The Pauli operator 

In this section we will study the Pauli operator corresponding to the magnetic 
field (J23J) obtained as the square of a self-adjoint Dirac operator from the previous 
section. 

Definition 3.1. We define the Pauli operator <$'' as (T) h ) 2 where T) h is a self- 
adjoint Dirac operator defined in Definition 12. II This means that 

V(<$ h ) ={^e {L 2 f : J>V € B(D' 1 )}, 
Vp h iP = (d h ) 2 ^, for V G T>(q$ h ). 
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Let us again introduce the boundary value linear functionals acting on 2)(<p' 1 ), 
but this time for all singular points A. For j = 1, . . . , n, let 



i c2tt 



c* Q ,(V±) = lim rf — / ^(rje^Wi, 



r,- Q;f c_ a .(r/'±) 



,.. ■. _ u_ !-<** J_ 



2;t 



c^-iGk) = Bmr""'- / V±(^)e^c% and 

cf_ Q ,(^ ± )=^un o r^- 1 (^£% ± (r J -e^Oe^^-r^- 1 C ±._ 1 (^ ± ) 

Since there are more self-adjoint Pauli extensions than Dirac extensions correspond- 
ing to our singular magnetic field, it is clear that not all Pauli operators can be 
obtained as the square of a self-adjoint Dirac operator. 

Proposition 3.2. For an arbitrary self-adjoint Pauli extension ^p, it is the square 
of some self-adjoint Dirac extension 1) h if and only if the following equations are 
satisfied for all ip G 2)(^J) 

- ,, \ =icoi{Tj/2) J , (3.1) 

-zcot(r i /2)-— i-r , (3.2) 



c+ a] (ij + ) = 0, and (3.3) 

%_i(V'-) = 0. (3.4) 

Proof. Let T) h be a given self-adjoint Dirac extension, and let ip belong to T)(T) h ). 
Then for some constants /ij we have 

i> ~ — 1 , for i = 1, ...,n 

2 V (l-e'^^Jrp+O^^) ; 

We want to find the next term in the asymptotical expansion for ip such that 

'(1 + e iT i)cr{l - a 3 )rf^ l e- ie ^ 



2 V (l-e^Ha.yp J 

for some constants Vj . This means that ip must have the asymptotics 

/i£(l + e -,)a(l - o^r^-'e"^ - ^(1 - e^')cr(a j - l)r) _a 'e-*' + 0(r?- at ') N 

V ?(l-e^M aj )rp'4(l + e^>KK+0(r; + °') , 

as tends to z,-, for j = 1, . . . , n. From this it follows that 



2 



c+_ 1 (^ + ) = ^(l + e^Ml-a J ), (3.5) 



4- ai W+) = -^(1 - e iT >(«j - 1), (3-6) 

c:a j W-) = f(l-e iTj Ma j ), and (3.7) 

c-ty_) = -^(l + e iT 'M-ai). (3-8) 
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Moreover 

c± ai (lM = 0, and (3.9) 

c^-iW-) = (3.10) 

since no such singular functions belong to CD(2)' 1 ). The coefficients 

c+ . ("0+ )= arbitrary, and (3.11) 

c i-u •(V'-) = arbitrary (3-12) 

since such terms disappear (near the singular point) when applying 1) h . D 

Remark. The definition of ( Jp h can be written as 

DpP") - {^ e (L 2 (M 2 )) 2 : 2 ^ g (i 2 (K 2 )) 2 ; 

equations l|3.1 |l -l|3.4 |l hold for all Vf- 

We see also that r D(Vp h ) is exactly the subset of T)(T) h ) for which also the condi- 
tion JH2J holds. D 

3.1. Spin-flip invariance and Zero-modes. 

Proposition 3.3. The only self-adjoint Pauli extensions fp 1 "-' 1 = (S) 1 "'^) 2 that are 
spin- flip invariant under the transform V are these where for all j = 1, . . . , n we 
have Tj = 7r/2 or Tj = 3n/2. 

Proof. The proof is the same as for the Dirac operators, see Proposition 12. 91 □ 

Theorem 3.4. If tj = t, j = l,...,n then the dimension of the kernel of < ^ h is 
given by 

f{|n-$|} ifr = 0, 
dimker^' 1 = < {|$|} «/r = 7T, 

I otherwise. 

Proof. This is a direct consequence of Theorem 12 . 1 21 since ker^ = kerS)' 1 . □ 

3.2. Discussion. Let us compare the different self-adjoint Pauli operators from J2] 
(which we will denote by *Pev) and d] (which we will denote by ^3 max ) with the 
ones obtained above as the square of a self-adjoint Dirac operator. It is easier to 
do this comparison if we have the same AB flux normalization for all operators. 
Thus, we let all AB intensities ctj belong to the interval (0, 1). In the case of the 
Pauli operator ^Pev, where the AB intensities were normalized to [—1/2, 1/2), we 
have to do a gauge transformation if there are intensities Oj belonging to [—1/2, 0). 
This is not a problem, since ^Pev is gauge invariant. 

In Table^we see a comparison of the boundary conditions of the Pauli operators 
obtained above that are the square of a Dirac operator and the Maximal and EV 
Pauli operators (see ^Hl El)- We see that ^p max is not the square of a Dirac 
operator. However, if we let 

(tt, if 0<a,<l/2 

r, ■ = < , , ? = 1, . . . ,n, 

J [0, if 1/2 < a, ; < 1 

and r = (n, . . . ,t„), then *Pev is the square of the self-adjoint Dirac operator 
corresponding to t. Note that it is possible to have different physical situations at 
the singular points A. Indeed, if not all intensities a.j belong to either (0, 1/2) or 
[1/2, 1) then this is the case. 
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Table 1. The boundary value conditions for the squared Dirac 
operators compared with the ones for the Maximal and EV Pauli 
operators. 





y$ h = {D h ) 2 


'Pniax 


*Pev 






3 

o± a . 


oo 


OO 


00 






c+ 


-gfe^fa/2) 





(oo, 


if 
if 


< a.j < 1/2 
1/2 < a, < 1 


c~ 

3 
a 3 


^^cot(r,/2) 





I ' 

[oo, 


if 
if 


< a, < 1/2 
1/2 < otj < 1 




oo 


00 


00 







Remark. If the AB intensities in J2| would have been normalized to (0, 1) instead 
of [—1/2, 1/2), then the operator *Pev would have become the square of the Dirac 
operator where tj — 7r for all j = 1, . . . , n. If the AB intensities would have been 
normalized to (—1,0) then ^Pev would have been the square of the Dirac operator 
where Tj = for all j = 1, . . . ,n. □ 

Among the Pauli operators studied in this article, the ones for r = tt/2 (which 
is (anti)-unitarily equivalent to the one for r = 37r/2), r = and r = 7r seems 
to be the most interesting ones. For r = 7r/2 we get a very symmetric domain 
of the operator, which implies that the operator is spin-flip invariant. Lacking 
zero-modes, it does not satisfy the original Aharonov-Casher formula, but it can be 
approximated component- wise according to Table ^ an d the result in [5] . See the 
end of 2U for a discussion of this. 

The Pauli operators corresponding to r = and t — n have very asymmetric 
domains. Only one of the components contain singular terms at the points A. This 
lack of symmetry implies that these extensions are not spin-flip invariant. On the 
other hand, the Pauli operator corresponding to r = n does satisfy the original 
Aharonov-Casher formula and there is no doubt that both of these Pauli operators 
can be approximated as in [8j, even as Pauli Hamiltonians. 

The Maximal Pauli operator studied ^H] is spin-flip invariant and has zero- 
modes, even more than is present in the original Aharonov-Casher formula. It 
can be approximated component-wise as in jHj. However, it is not the square of a 
self- adjoint Dirac operator. 

It is still not clear which Pauli extension that describes the physics in the best 
way. 
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